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Abstract 

We study a model initially proposed to describe a mixture of 
(CO) 1 _ a ,(N2) ;E adsorbed on exfoliated graphite. The approach used 
here is that of mean field theory. The Mean Field equations and the 
Helmholtz Free Energy are found. Phase diagrams are calculated too, 
and it is possible to find an analytic expression for the second order 
phase transition line. 

1 Introduction 

To calculate the fundamental equation of a many particles system, it 
is necessary first to determinate the allowed energy levels and make 
the sum of the partition function. Except for a small class of systems, 
that sum cannot be calculated exactly. One of the solutions for that 
problem is to look for approximated solutions, in our case we use the so 
called Mean Field Theory. The importance of this kind of calculation 
is that it gives us qualitative information about the critical behaviour 
of the system. 

Here we show results for a model of solid mixtures of molecules. We 
assume that one of the molecules have dipole and quadrupole moments 
(CO molecules) while the other one has only a quadrupole moment 
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(N2). This simple model is in good agreement with experimental re- 
sults [4]. The interactions between the CO molecules are assumed to 
be antiferroelectric, but between these and the N 2 molecules it is pro- 
posed a pseudo dipole-quadrupole coupling. The N2 molecules act like 
the random fields of the Random Field Ising Model. It is supposed 
that a crystal of polar molecules diluted with molecules which have 
quadrupole moments, can have a behaviour qualitatively similar to 
that predicted in random field models [5]. 

There are Monte Carlo results for this model [4, 6] which agree 
very well with experiments. Despite the fact that there are simula- 
tions and real experiments for this model, we don't have analytical 
results, exact or not, for this model yet. In order to improve our 
understanding about this kind of system, it is interesting to search 
for such results. The simplest way to find approximated results is to 
calculate a Mean Field Approximation for the model, what we do in 
the next sections. The paper is organized as follows. In section 2 we 
present the Dipole-Quadrupole model. In section 3 we calculate the 
mean field approximation using the BogoHubov Variational Theorem 
[1] obtaining the mean field equations and the Helmholtz Free En- 
ergy of the system. In section 4, we find analytic expressions for the 
transition lines and these lines are showed for several values of model 
parameters. 

2 Dipole-Quadrupole model 

This model was proposed to describe a mixture of (CO) 1 _ x (N2) a . ad- 
sorbed on exfoliated graphite [4]. We associate for each site of the 
lattice a spin Si = —1,0,1. Si = ±1 represents the electrical dipole 
orientations if the site is occupied by a CO molecule and Si = if 
the site is occupied by a N2 molecule. The interaction between the 
molecules of CO is antiferromagnetic and the interaction for CO - 
N2 pairs is described like a pseudo dipole-quadrupole one [4, 6]. The 
proposed Hamiltonian is 

n = - J E s * s i -JiY,^ 1 - s l )> w 

(ij) 

where J < represents the antiferromagnetic coupling, J\ describes 
the dipole-quadrupole interaction and Si = 0, ±1. The symbol (• • • ) 
denotes first neighbours summation. We rewrite the Hamiltonian in 
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the form 



H = - J ^ SiSjCiCj - Ji y~] SiCj(l - 9), 



(2) 



where Sj = ±1 and q = 1 if the site is occupied by a CO molecule 
and Ci = otherwise. The probability distribution for the occupation 
variables is 



P(a) =pS(c i -i) + (i-p)S(c i ), 



(3) 



with < p < 1. The presence of a antiferromagnetic interaction 
suggests that the system can spontaneously be subdivided in two sub- 
systems, represented by two sub-lattices (figure 1). 
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Figure 1: Sub-lattices A and B in a bidimensional lattice. 



3 Mean field calculations 

We propose an unperturbed Hamiltonian Hq of the form 
Ho = -( J\ + tja) ^2 SiCi ~ ( Jl + ^b) ^ SiCi > 

ieA i&B 



(4) 
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where ija and r\B are variational parameters and A and B indicate the 
inter-penetrating sub-lattices, each one with N/2 sites. We choose H$ 
to be, at the same time, soluble and similar in some manner to the 
original Hamiltonian. 

The Bogoliubov Inequality reads 



(5) 



(F) < (F ) + ((H-H ) t ), 



where Fq is the Helmholtz Free Energy calculated for the unperturbed 
system and F indicates the complete Helmholtz Free Energy. The 
symbol (• • • ) means a configurational average, while (•••)* indicates 
the thermal average weighted by the factors exp {—[3Ei) where Ei is 
the eigenvalue of energy with indice i. The partition function of the 
soluble Hamiltonian is 



z =TrK^}=nn 



X 



2 cosh y3ci(J\ + r)A 
2 cosh I (3cj{J\ + r\B 



(6) 



and the mean Free Energy calculated with that partition function 
becomes 



(Fo) 



N 
'2P 



In ^2cosh/3cj(Ji + tja)^ 
+ (2cosh/3Q(Ji +vb)) 



(7) 



the second part of inequality (5) is written 



(h - H Q )t) ={jJ2J2( s ^ c ^ + J i E T,& c * c ^ + 

ieAjeB ieAjeB 

r]A^2{SiCi) t + VB^2{SjCj) t y 
ieA jeB 



(8) 



If we adjust the variational parameters t\a and t]b using the stationary 
condition for the Free Energy and define the magnetizations 



mA = ( Cj tanh 



dp(H + f]A) 



(9) 
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and 



rris = ( Ci tanh 



dp(H + riB) 



(10) 



it is possible to write the Bogoliubov Inequality like 

In ^2 cosh (PaiJi + Va)))) + (in (2 cosh {f3c i {J l + rj B )) 
Ci)mA + (ci)m B 



v 1 ~ 2/3 



Nz Nz 
+ J—m A m B + J\— - 
2 4 



TV JV 
+ J/iymi + r] B —m B . 



(11) 



This expression must be minimized in order to estimate the Free En- 
ergy. Differentiating the equation and imposing the stationary condi- 
tion we find the variational parameters 



and 



rj A = -Jzm B - Ji-(ci) 



tie = -Jzm A - Jl-(Cj) 



(12) 
(13) 



then 



In ^2cosh (/3cj(Ji - -^-(q) - Jzm B ) 



+ (in (2 cosh ((3ci(Ji - -y-(cj) ~ Jzm A ) 



(14) 



— -Jzm A mB- 

Again, the antiferromagnetic interaction suggests we should define 

m A + m B 



M = 

2 

where M is the total magnetization and 

m A - m B 

m s = . 

2 

It is also convenient to define 



(15) 



(16) 



(17) 
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and 



Ji 
\J\z 



(l -§<«>), 



(18) 



we find 



(F) 



1 



+ 



In (^2 cosh (j(h - M + m,))) 
In ^2cosh(y(/i- M-m a ))) 



(19) 
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(15) and (16) define the mean field equations for this model 



M = - 
2 



and 



m., 



Cj tanh ^y (/i — M + m s )^ ^ + tanh ^y (/i 



Cj tanh ( y (/i — M + m s 



(20) 



Ci tanh ( y (h — M — m s 



(21) 



These coupled equations can be solved for M and m s using a iterative 
scheme. 



4 Critical lines 

The existence of a critical line, critical points and the determination of 
transition order are made by means of a Landau expansion [2] . Then 
we write the magnetization in the form 

M = M + m, (22) 

where Mq is a paramagnetic solution which is a solution of 

Mo = (citanh(^(/i-Mo))). (23) 



6 



We define the quantities m\ = m — m s and 7712 = m + m s to expand 
the right hand side of equation (22) near a paramagnetic soiution. 
Then 



M=Y,Mm 1 l + m 1 i), 

n=0 

eliminating the term corresponding to n = we find 

00 

n=l 

it is useful to expand the staggered magnetization too 

00 

m s = ^A n {ml-m%). 



n=l 



We determine the coefficients A n in these expressions 

Li = 4(c,(i-r a )), 



A, 



•>2 • -'/») )■ 



where 

T k = tanh fc 

It is possible to calculate m in function of m 



\{h-M Q ) 



(24) 



(25) 



M= 3^( C;3 ( 1_4T2+3T4 ))' 

A 4 = ^ i ^(3r 5 -5T 3 + 2r 1 )^ 



As = 15V ( c * 5 ( 15Te " 3m + lm " 2 ) ) ' (31) 



(26) 

(27) 
(28) 
(29) 
(30) 



(32) 



m = B\m 2 s + B2m s + B^m® + 



(33) 
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if we substitute this expansion in equation (25) and equal terms of 
same power in m s we find the expressions for the coefficients Bi in 
terms of A- 



B, 



A, 



l-^i ' 



(34) 



Bo 



A\ 



i-A 1 \(i-A 1 y 



+ Ai + 3 



A 2 A 3 
I -A, 



(35) 



Bo 



1-A, 



2Al 



A 3 



(1-^)2 1 (1-^)2 



AUi 



A 4 + 3 



A 2 A 3 
l-Ai 



i2^3 SA\M 5^ 2 ^ 5 

+ (l-^i) 3 + (1-A 1 )*1-A 1 
Here, substituting (33) in (26) we arrive at expression 

am s + bml + cm 5 s + • • • = 0, 
and finally at the coefficients 

a = -(l + A 1 ), 



A\ 



1-A, 



+ A 3 \, 



(36) 



(37) 

(38) 
(39) 



A 2 



A 3 



1-A 1 \(l-A 1 y 



+ A 4 + 3 



A 2 A 3 
I -A, 



+ 3 



AlA 



nSi 3 AoAd 

1 +4^-^ + ,45 



(l-^i) 



1-A 1 
(40) 



A Second order transition is determined by the condition a = and 
b > 0. The existence of a tricritical point will be determined by 
o = 6 = with c> [3]. 

With these conditions found, we calculate configurational mean 
values as indicated for the distribution eq. (3) to find that the critical 
line condition is equivalent to 



1 



l-tanh^(^) 



0, 



(41) 
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which, together with the configurational mean value of expression (23) 
gives us 



Mo 
P 



P 



1/2 



(42) 



from which we note a limitation on reduced temperature represented 
by t < p. Now the critical line equation is 



h = t tanh 



r t~ 


1/2 


r t- 


1/2 


l - - 


+ P 


l - - 




L p. 


p- 





(43) 



and, finally 

A _ 



;i-f) 



t tanh 



P 



(44) 



We should note here that, for a particular lattice of coordination num- 
ber z, there exist a critical value for the occupation p where there is 
a divergence. This point is located at 



Pc 



2 

z 



(45) 



The condition c > for existence of a tricritical point 



2 4 fh- M \ (2p p 1\ 2 fh- M \ p p 

(46) 

Is not satisfied for this model. 

In figures (2) and (3) we have two phase diagrams (one for p < 2/z 
and one for p > 2/z) for the model in the case z = 4 indicating two 
possible phases: Antiferromagnetic (AF) and Paramagnetic (P). In 
figures (4) and (5) we compare the phase diagrams for different values 
of occupation p. 
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Figure 4: Comparison for three values of the occupation rate which are less 
than the critical occupation. 




Figure 5: Comparison for three values of the occupation rate which are bigger 
than the critical occupation. 



11 



5 Conclusions 



Preliminary mean field results for the model Hamiltonian (2) was pre- 
sented in this paper. This model describes dipolar systems diluted 
with quadrupolar molecules. It was possible to determine, using the 
Bogoliubov Inequality, the mean field equations for the model and the 
Helmholtz Free Energy. The last one was calculated both, in terms 
of the magnetizations of the sub-lattices m a and mi, and in terms of 
total magnetization M and staggered magnetization m s . We have also 
calculated the analytical form for the second order jjr x t critical line. 
We found a critical value of occupation probability p c in which there 
is a divergence in expression (44). The critical lines are remarkably 
different in the regimes of p < p c and p > p c . Another interesting fact 
is that this model does not have any tricritical or multicritical points, 
at least in the mean field approximation. 

We hope this study is able to stimulate theoretical physicists to 
search for more theoretical results about this kind of system. It is 
specially interesting to learn about the possible connections between 
these and random field models. 
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